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The heavy quark-antiquark potential is accessible in perturbative QCD and in lattice simu- 
lations. The perturbative short-distance part of the potential is contructed via a restricted 
Fourier transform, covering the momentum region where perturbative QCD is applicable. 
We show that for the leading order static term as well as for the mass dependent corrections, 
the perturbative part can be matched at intermediate distances with results from lattice QCD. 
From these matched potentials, quarkonium spectra with a single free parameter (the heavy 
quark mass) are derived and compared with empirical spectra. Furthermore, charm and bot- 
tom quark masses are deduced. 



1 The static potential 

The potential between two heavy quarks is a prime subject of interest since the early days 
of QCD. Nowadays it is defined in a non-relativistic effective theory framework. While the 
long distance part can be studied in lattice QCD simulations, perturbation theory should 
be expected to work at short distances. The potential can be organized in a power series of 
the inverse quark mass m: 

(1) V = + -- — \ — h 

The leading term V^^^ represents the static potential. It has the following form at two-loop 
order in momentum space: 

(2) 17(«)(|,-|) = + + {"-WP^^- 
where q is the three-momentum transfer. The constants fli and fl2 are [1-3]: 

(3) fli = 31/3- 10/9 n/, 

(4) ai = 456.749 - 66.3542 rif + 1.23457 nj, 
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Figure 1 : Static QCD potential (with rif = 3) from the restricted numerical Fourier trans- 
form (6). Left: coordinate space potential at two-loop order for different values oi }ijr. 
The curves have been shifted by a constant to match at small r values. Right: potential 
matched at r = 0.14 fm to a potential from lattice QCD [7]. Taken from Ref. [8]. 



where is the number of light-quark flavors. Higher order terms have infrared contribu- 
tions and are not considered at this level. Expressing ( | (f" | ) in a power series expansion 
about Ks at a fixed scale }i leads to the standard definition of the r-space static potential. 



(5) y(°)(r) 



47r 



47r 



with gii{r) = ln(/^r) + 7E . It is well known that this potential suffers from renormalon 
ambiguities [4,5] and shows a badly convergent behavior [6]. 

We work in the following in the potential subtracted (PS) scheme proposed by Beneke [4] 
and define the static r-space potential. 



(6) 



\<l\>l'f 



(271)3 



where V'^^H\q\) is given in Eq. (2), but a, ( |(f | ) is understood without resorting to a power 
series expansion. The momentum space cutoff }if is introduced in order to exclude the 
uncontrolled low momentum region. The potential V^^\r,}ij^) is evaluated numerically 
using four-loop RGE running for the strong coupling ag- For distances r < 0.2 fm, this 
potential depends only marginally on /^^ as shown in the left plot of Fig. 1. The perturba- 
tive potential, valid at small distances, can be matched at intermediate distances to results 
from lattice QCD (see the rightmost plot in Fig. 1). For the matching point (dashed line) 
we choose r = 0.14 fm where both the perturbative and lattice potential are expected to be 
reliable. 
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Figure 2: The order 1/m potential with = 3, from the restricted Fourier transform. 
Left: perturbative potential for different cutoffs Right: perturbative potential matched 
at intermediate distances to a potential from lattice QCD. Taken from Ref. [8]. 



2 The order 1/m potential and quarkonium spectroscopy 

in Eq. (1) is the first mass dependent correction to the static potential. It is spin inde- 
pendent and the leading term reads in momentum space [9]: 

(7) t7W(|,n) = -^^^||^{l + 0(..)}. 

It can be transformed analogously as in Eq. (6) to r-space with a low momentum cutoff 

The dependence of V'^^'> on the cutoff scale is again very weak for distances r < 0.2 fm. At 
long distances y'^' (r) is known from lattice QCD [10, 11]. To fit the lattice data we use the 
form 

(8) vm^,~^^_c_^^,^^^L^^^ const, 

motivated in [12]. As shown in Fig. 2 matching with the perturbative potential at interme- 
diate distances is also possible at order 1/m. 

Using y'*^' and V^^'^ as input in the Schrodinger equation, we can examine bottomonium 
and charmonium spectra. The overall constant of the potential is the only free parameter. 
This single parameter is related to the heavy quark mass in the PS scheme and can be 
translated in a second step to the bottom and charm quark masses in the MS scheme (see [8] 
for details). Our findings for the masses are summarized in Table 1 and compared to the 
values listed by the Particle Data Group (PDG) [13]. Results for the bottomonium and 
charmonium spectra are shown in Fig. 3. In both cases we find that the IS states are the 
most strongly affected by 1/m-effects. An additional effective one-gluon exchange spin 
dependent term with af^ = 0.3 (+h.f. in Fig. 3) would improve our predictions. Of course. 
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MS masses [GeV] 





Static 


Static + 


PDG 2010 


Bottom quark 


4.20 ± 0.04 


4 1 8+0.05 


4 19+0.18 


Charm quark 


1.23 ± 0.04 


1 28+°-°^ 


-1 97+0.07 
^•^'-0.09 



Table 1 : Comparison of quark masses obtained in our approach (leading order plus or- 
der 1/m corrections) with the values listed by the Particle Data Group (PDG) [13]. 

this step is purely ad hoc and r\eeds to be substituted by the full potential of order l/rrp-, 
to be investigated in forthcoming work. 
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Figure 3: Bottomonium and charmonium spectrum in comparison with experiment. 
Static plus order 1 / m results are shown, with additional hyperfine effects (h.f.) added 
phenomenologically. Taken from Ref. [8]. 
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